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Abstract 

The next-next-to-leading order QCD corrections to the e + e~ annihilation into 
hadrons are considered. The stability of the predictions with respect to change 
of the renormalization scheme is discussed in detail for the case of five, four 
and three active quark flavors. The analysis is based on the recently pro- 
posed condition for selecting renormalization schemes according to the degree 
of cancellation that they introduce in the expression for the scheme invariant 
combination of the expansion coefficients. It is demonstrated that the scheme 
dependence ambiguity in the predictions obtained with the conventional ex- 
pansion is substantial, particularly at lower energies. It is shown however, 
that the stability of the predictions is greatly improved when QCD correc- 
tions are evaluated in a more precise way, by utilizing the contour integral 
representation and calculating numerically the contour integral. 
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1 Introduction 



In a series of recent papers 0-fl a method has been presented for a systematic 
analysis of the renormalization scheme (RS) ambiguities in the next-next-to-leading 
(NNLO) perturbative QCD predictions. It was emphasized in [Q, fl that besides 
giving predictions in some preferred renormalization scheme one should also inves- 
tigate the stability of the predictions when the parameters determining the scheme 
are changed in some acceptable range. The method discussed in [[!], Q involves a 
specific condition that allows one to eliminate from the analysis the renormalization 
schemes that give rise to unnaturally large expansion coefficients. The condition 
on the acceptable schemes is based on the existence in NNLO of the RS invariant 
combination of the expansion coefficients, which is characteristic for the considered 
physical quantity. The method of |], |2|, [3| has been applied to the QCD corrections 
to the Bjorken sum rule for the polarized structure functions |§ and to the QCD 
corrections to the total hadronic width of the tau lepton ^ [|. 

In this note we apply this method to the QCD correction to R e + e - ratio: 

o tot (e + e- -> hadrons) 

K e+e 1 T _ T _s ■ (1J 



which received considerable attention in recent years |E|-[27[|. We show that a 
straightforward application of the condition proposed in Jl], [|, |3| to the conventional 
perturbative expression for the QCD effects in the R e + e - ratio exhibits a rather 
strong RS dependence, even at high energies. Looking for improvement and moti- 
vated by the analysis of the corrections to the tau decay JT], |2], [| , we calculate 
the QCD correction to the R e + e - ratio by using the contour integral representation 
pn| and evaluating the contour integral numerically. In this way we resumm to 
all orders some of the so called 7r 2 corrections, which appear as a result of analytic 
continuation of the expression for the hadronic vacuum polarization function from 
spacelike to timelike momenta |32|, Such corrections constitute a dominant con- 
tribution in the NNLO. Using the improved expression we perform similar analysis 
as in the case of the conventional expansion. We find that the predictions obtained 
by numerical evaluation of the contour integral show extremely good stability with 
respect to change of the RS. 

The results reported here have been announced in || and briefly described in 
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2 S e + e - and the problem of renormalization scheme 
ambiguity 

Away from the thresholds, neglecting the effects of the quark masses and the elec- 
troweak corrections, the formula for R e + e - may be written in the form: 

R e+e -(s) = 3^2Q 2 f [l + 5 e+e -(s)}, (2) 
/ 

where Qf denotes the electric charge of the quark with the flavor / and S e + e - is the 
QCD correction. The renormalization group improved NNLO expression for S e + e - 
has the form: 

^? e _(s) = a(s) [l+r ia (s)+r 2 a 2 (s)}, (3) 

where a(/i 2 ) = g 2 (/i 2 )/(47r 2 ) is the coupling constant, satisfying the renormalization 
group equation: 

a— = —ba 2 (1 + Ci<2 + c 2 a 2 ). (4) 

(2) 

The perturbative result for o e + e - is usually expressed in the Modified Minimal Sub- 
traction (MS) renormalization scheme, i.e. using the MS renormalization convention 
with fi 2 = s. In the MS scheme we have ^5|, 0, g: 



rf s = 1.985707- 0.115295 n f , (5) 
rfs = 18.242692- 4.215847 n/ + 0.086207nj + 

+ rr 9 -(W2) 2 A (6) 

where the r^™ 9 term in r^f s represents the so called flavor singlet contribution: 

(EfQf) 2 ( 55 5 



which arises from the light-by-light scattering type of diagrams (£3 = 1.202056903). 
(It should be noted that the first calculation of the NNLO correction was erro- 
neous. The corrected result was published in 0. An independent evaluation was 
reported in 0.) For the coefficients in the renormalization group equation we have 
b = (33 - 2n/)/6, c x = (153 - 19n/)/(66 - 4n/) and §§: 



77139 - 15099 n f + 325 n 2 

cf = ; -■ 8 

2 288(33- 2 nf) V ; 

For convenience we collect in Table 1 the numerical values of the expansion coeffi- 
cients for various values of n/. 

Besides the MS scheme other choices of the RS are of course possible, such as for 



example the momentum subtraction schemes [37|. A change in the RS modifies the 
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r MS 
' 1 


r MS 
'2 


smg 
r 2 


Cl 


r MS 
c 2 


P2 


2 


1.75512 


-9.14055 


-0.08264 


1.98276 


5.77598 


-9.92498 


3 


1.63982 


-10.28394 


0.00000 


1.77778 


4.47106 


-11.41713 


4 


1.52453 


-11.68560 


-0.16527 


1.54000 


3.04764 


-13.30991 


5 


1.40923 


-12.80463 


-0.03756 


1.26087 


1.47479 


-15.09262 


6 


1.29394 


-14.27207 


-0.24791 


0.92857 


-0.29018 


-17.43803 



Table 1: Numerical values of the expansion coefficients r« for S e + e - } obtained with 
the MS renormalization convention and p? = s, for various numbers of quark flavors. 
The magnitude of the flavor singlet contribution r^ 719 is separately indicated. The 
values of the RS invariant p% are calculated according to Eq. (^). The numerical 
values of the coefficients q in the renormalization group equation are included for 
completeness. 



values of the expansion coefficients — the relevant formulas have been collected for 
example in 0. (The coefficients b and c\ are RS independent in the class of mass and 
gauge independent schemes.) The change in the expansion coefficients compensates 
for the finite renormalization of the coupling constant. Of course, in the given 
order of perturbation expansion this compensation may be only approximate, so 
that there is some numerical difference in the perturbative predictions in various 
schemes. This difference is formally of higher order in the coupling — it is 0(a 4 ) 
for the NNLO expression — but numerically the difference may be significant for 
comparison of theoretical predictions with the experimental data. There has been a 



lively discussion how to avoid this problem, both in the general case f3~8]|- ||4"I|| (for a 
summary of early contributions see |42]| ) and in the particular case of S e + e - pT|-l|20||. 
(Some of the early papers [p[| -[p!5| contain discussion of 8f+ e - with the erroneous 
value of the NNLO correction reported in M. Much of the initial interest in the 

(2) 

RS dependence of o + - came from the fact that this erroneous correction was very 
large.) It seems that one of the most interesting propositions is to choose the scheme 
according to the so called Principle of Minimal Sensitivity (PMS) |39| j. 

However, as was emphasized in [p], |^, |3|, besides calculating the predictions in 
some preferred renormalization scheme, it is also important to investigate the stabil- 
ity of the predictions with respect to reasonable variations in the scheme parameters. 
By calculating the variation in the predictions over some set of a priori acceptable 
schemes one obtains a quantitative estimate of reliability of the optimized predic- 
tions. A systematic method for analyzing the stability of predictions with respect to 
change of the renormalization scheme has been presented in [jl], |2|, [I| . This method is 
based on the existence of the RS invariant combination of the expansion coefficients 

p 2 = c 2 + r 2 - can - r 2 , (9) 
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which appears to be a natural RS independent characterization of the magnitude 
of the NNLO correction. (We adopt here the definition of the RS invariant used in 
|38| , f4Tf , which differs by a constant from the definition of Stevenson |39|]: p 2 tev 



p 2 — cf/4. The arguments in favor of Eq. have been given in ||.) The numerical 
values of this invariant in the case of S e + e -, for different values of rif, are collected 
in Table 1. 

The p2 invariant may be used to eliminate from the analysis the unnatural renor- 
malization schemes. This is done by introducing a function a 2 defined on the space 
of the expansion coefficients: 

^2(^1,^2,02) = |c 2 | + |r 2 | +ci|ri| + r\, (10) 

which measures the degree of cancellation in the expression for p 2 . An unnatural 
renormalization scheme, which artificially introduces large expansion coefficients, 
would be immediately distinguished by a value of a 2 which would be large compared 
to \p 2 \. The function cr 2 defines classes of equivalence of the perturbative approx- 
imants. If one has any preference for using a perturbative expression obtained in 
some optimal scheme, one should also take into account predictions obtained in the 
schemes which imply the same, or smaller, cancellations in the expression for p 2 , i.e. 
which have the same, or smaller, value of a 2 . In particular, for the PMS scheme 
we have a 2 w 2\p 2 \ Q. Therefore it appears that the set of schemes which generate 
approximants satisfying o 2 < 2\p 2 \ is a minimal set that has to be taken into account 
in the analysis of stability of the predictions with respect to change of the RS. More 
generally, it is useful to use a condition on the allowed schemes in the form: 

<72{ri,r 2 ,c 2 ) < I \p 2 \, (11) 

where / > 1 is some constant, which determines how strong cancellations in the 
expression for p 2 we want to allow. 

In this note we analyze the RS dependence of the NNLO predictions for 5 e + e -, 
using systematically the condition ([11]). As in the previous papers |], 0, [|, we 
use the ri and c 2 coefficients as the two independent parameters characterizing 
the freedom of choice of the approximants in the NNLO. To obtain the numerical 
value of the running coupling constant we use the implicit equation, which results 
from integrating the renormalization group equation (Q) with appropriate boundary 
condition [34] : 

= r p_ ri + $( a;C2 ) ; (12) 




where 

$(a, c 2 ) = ci In ( -M + - + ci ln(cia) + 0(a). (13) 
\2ci J a 

The explicit form of <l>(a, c 2 ) is given for example in P5| . The appearance of A^- and 
r MS j n ^ e expression (O) is a result of taking into account the so called Celmaster- 



Gonsalves relation [37] between the lambda parameters in different schemes. This 



relation is valid to all orders of perturbation expansion. 
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mm 
' 1 



max 
I 1 



The region of the scheme parameters satisfying Eq. (|TT| ) has simple analytic 
description. In the case p 2 < and |p 2 | > 2c\{l + 1)/(Z — l) 2 let us define: 

-J\p 2 \(l + l)/2, (14) 

[_ Cl + v / c 2 + 2 (/ + l)|p 2 |]/2, (15) 

cT n = -M(; + l)/2, (16) 

cT x = |P2|(i-l)A (17) 

4"* = cir™ n + c™ aa; - (18) 

For c 2 > the region of allowed parameters is bounded from above by the line 
joining the points (r^ ir \0), (r^ in , cf l ) , (0,c^ aa; ), (rf oa: , d? ax ), (rf^O). For c 2 < 
the region of allowed parameters is bounded from below by the lines: 

c 2 ( ri ) = r 2 + cf n for rf n < ri <0, (19) 
c 2 ( ri ) = r 2 + can + c™ in for < r x < r™ ax . (20) 

In the case p 2 < and |p 2 | < 2c 2 (/ + l)/(/ — l) 2 we should use instead: 

rT n = -\P2\(1-1)/(2 CI ), (21) 

c int = ^mm^2 + c mm_ ( 22 ) 

(23) 

For c 2 > the region of allowed parameters is then bounded from above by the line 
joining the points (r mm ,0), (0,c^' ax ), ( r max ^ c max ^ (r max ,0). For c 2 < the region 
of allowed parameters is bounded from below by the line joining the points (r™ in , 0) 
and (r™ 1 ™, c 2 n *), and the curves defined in the previous case. 

For p 2 > c 2 /4 the region of the scheme parameters satisfying the Eq. ( |11~D has 
been described in ||. 

3 Estimate of the RS ambiguities in the conven- 
tional expansion for S e + e - 

Let us first consider the case of five active quark flavors, which is most important 
for experimental determination of A^r. The same corrections gives also a dominant 
QCD contribution to the hadronic width of the Z° boson. For nj = 5 we have 
P2 = —15.09262. In Fig. 1 we show the contour plot of 5^ as a function of the 
parameters r\ and c 2 , for y^i/A^- = 75. We have indicated the region of parameters 
satisfyings the condition ( |TT| ) with / = 2. For comparison, we also indicate the region 
corresponding to / = 3. The PMS prediction is represented in Fig. 1 by a saddle 
point at T\ = —0.408 and c 2 = —23.154. We see that the PMS parameters are close 



to the approximate solution of the PMS equations [44 



r? MS = 0(a PMS ), c p 2 MS = \p 2 + 0(a PM5 ), (24) 
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and the PMS point lies indeed on the boundary of the I = 2 region, as expected. 
Comparing the values of 5^ obtained for the scheme parameters in the I = 2 
region we find for y^/A^^ = 75, that the minimal value is attained for r± = —4.76, 
C2 = 1.55 and the maximal value is attained for n = 3.52, C2 = 7.55. For the I = 3 
region we obtain the minimal value for n = —5.49, C2 = 8.17, and the maximal 
value for r% = 3.98, C2 = 15.09. In both cases the maximal and minimal values are 
attained at the boundary of the allowed region. Let us note, that the commonly 
used MS scheme lies within the I = 2 region. 

Performing similar contour plots in the range 40 < y/s/A^ < 200 we find, 

(2) 

that the scheme parameters, for which o + reaches extremal values in the I = 2, 3 
allowed regions, are practically independent of the y^/A™-. 

(2) 

In Fig. 2 we show how the maximal and minimal values of o + in the I = 2, 3 
allowed regions depend on \/s/Aj^. We also show the PMS prediction and the 
experimental constraint 8%? e -(y/s = 31.6 GeV) = 0.0527 ± 0.0050 |TJ. We find 
that with increasing y/s/Aj^ the scheme dependence is decreasing, as expected, 
although it remains substantial even for high energies. Let us take for example 
y/s/ Ajj^ = 162, corresponding to A|^ = 0.195 GeV — which is the value preferred 
by the Particle Data Group [HH] — and \fs = 31.6. In this case the scheme variation 

(2) 

of 5 + over the 1 = 2 region is 5% of the PMS prediction, and for the 1 = 3 region 
8%, compared with 9% and 16%, respectively, for y/s/ Aj^ = 75. However, when we 
decrease y/s/Aj^ below 75 the scheme dependence increases rapidly, and it becomes 



very large already for y/s/A^^ = 30. The scheme dependence appears to be quite 

(2) 

large in the range of values of o + e - relevant for fitting the experimental data. For 
example, the line representing the minimal values on the I = 2 region does not 
reach the central experimental value, which translates into a very large theoretical 
uncertainty in the fitted value of Ajjg. 

For rif = 4 we have pf = —13.30991. In Fig. 3 we show the contour plot of 
5g+ e _ as a function of the parameters r\ and C2, for ^Js/Ajj- = 30. Similarly as 
in the rif = 5 case we find that the PMS prediction is well represented by the 
approximate solution (^5). The variation over the I = 2 region is approximately 
11% of the PMS prediction. In Fig. 4 we show the variation in the predictions for 

(2) 

Og +e _ when the scheme parameters are changed over the I = 2 region, as a function 
of ^/s/A-^,. It is evident that for y^/A-^ smaller that 20, which is the range 
relevant for fitting the experimental data, this variation becomes very large. (Note 
that analysis of experimental data from several experiments gives [^5[ 5l+ p e -(y/s = 
9 GeV) = 0.073 ±0.024.) 

Finally, for nj = 3 we have = —11.41713. In Fig. 5 we show the contour plot 
of 8\2- as a function of the parameters r\ and C2, for y/s/Aj-f^ = 9. The variation 
of the predictions over the / = 2 region is approximately 28% of the PMS value. In 

(2) 

Fig. 6 we show the variation in the predictions for 5 + when the scheme parameters 
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are changed over the / = 2 region, as a function of y/s/A-^. We observe that the 



.(3) _ 



variation in the predictions starts to increase rapidly for y^/Ai^L smaller than 9. 

Let us summarize our analysis of the predictions for 5 e + e - obtained from the 
conventional expansion. We found that changing the renormalization scheme within 
a class of schemes which, according to our condition (|TTD , appear to be as good as 
the PMS scheme, we obtain rather large variation in the predictions. In some cases 
we may even speak about instability of the predictions with respect to change of 



the renormalization scheme. This is in contrast with the statement in ||16|| , that the 
conventional expansion for 5 e + e - is highly reliable. The conclusion found in |16] is 
based on the observation, that for 5 e + e - the MS prediction is very close to the PMS 
prediction. The fact that the MS prediction is very close to the PMS prediction is 
of course true — for example in the scale of Fig. 2 the MS and PMS curves would 
be difficult to distinguish. Similar situation occurs for other values of rif. It is 
clear however, that there is no theoretical or phenomenological motivation to use 
the MS-PMS difference as a measure of reliability of the perturbation expansion for 
any physical quantity. The fact that the MS prediction for S e + e - is close to the PMS 
prediction is simply a coincidence, without deeper significance for such problems 
as reliability of the predictions and good or bad convergence of the perturbation 
expansion. 

It is interesting to note that for very low energies the PMS predictions display 
the infrared fixed point type of behavior ||18|. However, this type of behavior, which 



in fact does not manifest itself in the rif = 3 predictions until y^s/ Aj^ m 2.5, is ac- 
companied by a rapidly increasing RS dependence. It seems therefore unreasonable 
to put too much faith in the PMS prediction when even a very small change of the 
scheme parameters dramatically modifies the result. These remarks apply as well 
to the case rif = 2. 



4 Analysis of the tt 2 terms in 5 e + e - 

The strong RS dependence described above is somewhat surprising. It may seem 
understandable that the perturbation expansion is not reliable in the energy range 
appropriate for example for the rif = 3 regime. However, one would expect that 
y^s/A^^ of order 75 is large enough for the perturbation series to be very well 
behaved. The origin of the strong scheme dependence may be traced back to the 
fact that the NNLO correction is relatively large, which is reflected by large value of 
the RS invariant However, a major contribution to the NNLO correction comes 
from the term which appears in the process of analytic continuation of perturbative 
expression from spacelike to timelike momenta. To see clearly the significance of 
such contributions, and to show how one may treat them in an improved way, it is 
convenient to use the so called Adler function pTf] : 

D(g 2 ) = -127rV-^n(g 2 ). (25) 
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where II(g 2 ) is the transverse part of quark electromagnetic current correlator IP^(g): 

n^(g) = {-<rtf + <f<nntf), (26) 



W u (q) 



d 4 xe iqx <0\T(J fl (x)r(0y)\0> 



Neglecting the quark mass effects and electroweak corrections we may write: 

f 



(27) 



(28) 



where 8d{— q 2 ) denotes the QCD correction. The Adler function is RS invariant in 
the formal sense, i.e. it may be considered to be a physical quantity, despite the 
fact that it cannot be directly measured in the experiment. In particular, Sd{— q 2 ) 
is renormalization group invariant, in contrast to II(g 2 ), which does not even satisfy 
a homogenous renormalization group equation [fjSfl . The Adler function is directly 
calculable in the perturbation expansion for spacelike momenta. To express the 
R e + e - ratio by the Adler function one inverts the relation ([2~5|): 



n(g 2 ) - i% 2 ) 



127T 2 



da 



D(a) 



a 



(29) 



where q 2 is some reference spacelike momentum, and one utilizes the relation: 

6n 



R e+e . ( s ) = 12 7r Imn(s + ie) 



[U(s + ie) -Il(s-ie)] 



(30) 



In this way one obtains R e + e -(s) as a contour integral in the complex momentum 
plane, with the Adler function under the integral pi BE: 



(31) 



1 f , D(a) 
2-ni Jc a 



where the contour C runs clockwise from a = s — ie to a = below the real positive 
axis, around o = 0, and to a = s + ie above the real positive axis. The integration 
contour may be of course arbitrarily deformed in the domain of analyticity of the 
Adler function. A convenient choice is q 2 = —sexp(—i6) with 9 G [— tt, tt]. For this 
choice of the contour we obtain the following simple relation between 5 e + e - (s) and 
5 D (-q 2 ): 

1 



8 e 



2tt 



d6 



Sd(-(t) 



5 exp(— i6) 



(32) 



The conventional expression for 5 e + e -(s) may be recovered from this formula by 
inserting under the contour integral an expansion of 5r>(— q 2 ) in terms of a(s): 



4 2) ("<7 2 ) = a(s) {l+[h~ (6/2) H-q 2 /s)} a(s) + 



+ 



h ~ (ci + 2fi)(6/2) H-q 2 /s) + (6/2) 2 (ln(-g 2 / S )) 2 ] a 2 (s)} , (33) 
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n f 




2 


9.28877 


3 


5.23783 


4 


0.96903 


5 


-3.00693 


6 


-7.36281 



Table 2: Numerical values of the RS invariant p® characteristic for the QCD 
correction to the Adler function. 



where fj denote the coefficients for expansion of Sd{—Q 2 ) in terms of a(—q 2 ). Eval- 
uating the trivial contour integrals involving powers of m(— cr/s), we obtain the 
expression @ with: 



1 bn" 

ri = ri, r 2 = 7*2 - - — 



(34) 



This implies pf = p^ — (^ 7r /2) 2 /3. In Table 2 we list the values of p^ for various 
values of n/. 

Numerically the contribution of the n 2 term is very large — for example for 
n f = 5 we have pf - pf = -12.08570. 



Contributions proportional to n appear also in higher orders. We have |25| : 

2 



r 3 



r 3 



bn 



(2f 2 + I Cl f 1 + ic? + c 2 ) 



1 /6vr 



(35) 
(36) 



The result for rs may be found in [22 



1 



r 5 



27 



7 



7*5 - - 10r 3 + — cif 2 + 4qri + -Cic 2 + 8c 2 f x + -c 3 — + 



&7T 
T 



1 77 

+ - 07*1 H Ci 

5 V 12 



bix 



(37) 



(The difference between r» and in higher orders was studied in |25], |33[.) Note 
that the 7r 2 corrections to 7*3 and are fully determined by the NNLO expression 
for 5d(— q 2 ). Taking into account that we have the following expressions for the 



higher order RS invariant combinations of the expansion coefficients [41 



p 3 = c 3 + 2r 3 - 4r 2 ri - 2rxp 2 - r^i + 2r 1; 
Pa = C4 + 3r 4 — 6r 3 ri — 4r^ — 3rip 3 — 4r^ — 
-(r 2 + 2r^)p 2 + \\r 2 r\ + Ci(r 3 - 3r 2 r x 



(38) 
(39) 
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we obtain: 




p? = If I , i m 

- ^-^? + ^(y) 2 4(t)*- <«> 

The 7r 2 terms are quite sizeable numerically. For example for rif = 5 we have: 

p£-p£ = -76.1924, p$ - p% = 211.025. (42) 

It is evident that the terms arising from the analytic continuation would make 
a significant contribution to the RS invariants in any order of the perturbation 
expansion. 

Returning to the evaluation of S e + e -(s), we note that the procedure used to 
obtain the conventional result treats the q 2 dependence of 5p in the complex energy 
plane in a rather crude way. A straightforward way to improve this evaluation is 
to use under the contour integral the renormalization group improved expression 
for 6d(— cr), analytically continued from the real negative a to the whole complex 
energy plane cut along the real positive axis. In other words, one should take into 
account the renormalization group evolution of a(—a) in the complex energy plane, 
avoiding the expansion of a(—a) in terms of a(s). In this way one makes maximal 
use of the renormalization group invariance property of the Adler function. Of 
course the integral may be now done only numerically, and the resulting expression 
for 5 e + e -(s) is no longer a polynomial in a(s), despite the fact that only the NNLO 
expression for the Adler function is used. It is easy to convince onself that the 
procedure outlined above is equivalent to the resummation — to all orders — of the 
7r 2 terms that contain powers of b, c\ and/or c%. (The summation of the leading 
terms proportional to (&7r/2) 2 was discussed in |33fl .) 

The improved approach based on the contour integral has been implemented with 
success in the case of the QCD corrections to the tau lepton decay [^R £|], where a 
similar problem of strong renormalization scheme dependence appears. It was found 
that using the contour integral representation and evaluating the contour integral 
numerically one obtains considerable improvement in the stability of predictions 



with respect to change of RS [£9, BJ. It is therefore of great interest to see whether 



one may improve in this way the predictions for S e 



5 Improved evaluation of S e + e - 

In this section we perform an analysis similar to that in the Section 3, using now the 
improved predictions for S e + e -, obtained by evaluating numerically the contour inte- 
gral in Eq. fl3"2"|). Similarly as in the case of the conventional perturbation expansion, 
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(2) 

we begin with the rif = 5 case. To show, how the improved evaluation of o y + af- 

(2) 

fects its RS dependence, we compare the plots of 0g+ e - as a function of ri, for several 
values of C2, with y/s/Aj^g = 75, obtained with the conventional NNLO expression 
(Fig. 7) and with the numerical evaluation of the contour integral (Fig. 8). We see 
that the predictions obtained by the numerical evaluation of the contour integral 
have much smaller RS dependence. In Fig. 8 we have also indicated the predictions 
obtained with the conventional expansion supplemented by the 0(a 4 ) and 0(a 5 ) 
terms given by Eq. fl3"5]) and Eq. fl36|). We see that this type of simple improve- 
ment of the conventional expansion reproduces quite well the results obtained from 
exact contour integral, except for large negative r\. (Inclusion of only the 0(a 4 ) 
term does not give good approximation. Inclusion of the 0(a e ) correction given by 
Eq. (^), which is of course only partially known at present, slightly improves the 
approximation for positive r\.) 

In Fig. 9 we show the contour plot of ^+ e - obtained from the expression (|32|) 
for y^/A^- =75. In Fig. 9 we also show the relevant regions of the scheme 
parameters satisfying the condition (|TT|) with I = 2,3. These regions are calculated 
assuming p 2 = pf because the basic object in the improved approach is 5$. For 
rif = 5 we have = —3.00693, which is much smaller in absolute value than p R . 
Consequently, the region of the allowed scheme parameters is much smaller than in 
the analysis of the conventional NNLO approximant. The improved predictions for 
5 e + e - have a saddle point type of behavior as a function of T\ and c^, where the saddle 
point represents the PMS prediction. However, the location of the saddle point is 
completely different than in the case of conventional expansion. (The location of 
the saddle point for the improved expression is no longer a solution of the set of 
the PMS equations given in [p9[] , because the improved approximant ( p2[ ) is not 
a polynomial in the running coupling constant.) It is interesting that the PMS 
point for the improved expression lies very close to the point ri = and c 2 = 
1.5pf = —4.51, which corresponds to the approximate value of the PMS parameters 

(2) 

if 5 D is optimized for spacelike momenta. Let us note that the MS scheme lies 
outside the / = 2 region in this case. However, the MS prediction in the improved 
approach is very close to the improved PMS prediction: we have 0.05279 and 0.05275 
respectively. 

The variation of the predictions over the I = 2 region is 0.3% of the PMS predic- 
tion, and variation over the I = 3 region is 0.5% of the PMS prediction. Even if we 
take variation over the region corresponding to / = 10 we obtain only 2.5% change 
in the predictions. We see that the improved prediction for 5 e + e - shows wonderful 
stability with respect to change of the RS. From Fig. 7 and Fig. 8 it is also clear, 
that the difference between NNLO and NLO PMS predictions is much smaller in 
the case of the improved prediction — 0.9% of the NNLO result for y^s/A^^ = 75 
- than in the case of the conventional expansion — 4.7% of the NNLO result. 
We conclude therefore that the theoretical ambiguities involved in the evaluation 
of 8 e + are in fact very small, provided that the analytic continuation effects are 
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^opt,NNLO 


^opt,NLO 


25 


0.06799 


0.06888 


50 


0.05753 


0.05811 


75 


0.05275 


0.05320 


100 


0.04981 


0.05019 


200 


0.04389 


0.04415 


500 


0.03791 


0.03809 



Table 3: Numerical values of the optimized predictions for S e + e -, obtained from 
the contour integral expression ([32]) for rif = 5. The PMS parameters are well 
approximated by r 1 = 0, c 2 = 1.5pf (NNLO) and r 1 = -0.59 (NLO). 

treated with appropriate care. For completeness, we give in Table 3 the NNLO and 
NLO PMS predictions in the improved approach for several values of y^s/Aj^. 

In the case of n/ = 5 predictions it is interesting how the improved eval- 
uation affects the fit to experimental data. Using the experimental constraint 
5j* e -Wa = 31.6 GeV) = 0.0527 ± 0.0050 |5[ and the improved PMS prediction 
we find A^^ = 0.419 ± 0.194 GeV, which is equivalent in the three loop approxi- 
mation to af 5 (M|) = 0.1319 ± 0.0100. For comparison, using the conventional 
expansion in the MS scheme we obtain the central value of A|^ = 0.399 GeV 
(a^ (M§) = 0.1308), while with the PMS prescription in the conventional ex- 
pansion we get AgL = 0.410 GeV (af^(M|) = 0.1314). We see therefore that 

(2) 

improvement in the evaluation of o e + has small effect on the fitted values of the 
A|^ parameter. 

For rif = 4 we have = 0.96903, i.e. the effect of tt 2 corrections is even larger 
than in the rif = 5 case. The rif = 4 case is in all respects similar to the rif = 5 case, 
except for the fact that the reduction in RS dependence seems to be even stronger. In 
Fig. 10 and Fig. 11 we compare the plots of <5g+ e - as a function of r±, for several values 

of c 2 , with y^/A^- = 30, obtained with the conventional NNLO expression (Fig. 10) 
and with the numerical evaluation of the contour integral (Fig. 11). In Fig. 11 we 
also show the predictions obtained with the conventional expansion supplemented 
by the 0(a 4 ) and 0(a 5 ) terms given by Eq. fl3"o]) and Eq. fl3"6"|). (Inclusion of the 
0(a 6 ) correction (|37|) does not improve the approximation.) In Fig. 12 we show the 
contour plot of the improved prediction for 5 e + e - obtained for y/s/Aj^ = 30. It 
is interesting that variation of the predictions over the I = 2 region is extremely 
small, of the order of 0.03% (!) of the PMS prediction. The improved prediction for 
v/s/A^r = 30 in the MS scheme is 0.05902, quite close to the improved PMS result 
0.05907. The differences with the results obtained in the conventional approach 
again are not very big — using the conventional expansion we have 0.06025 in the 
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front, N N LO 
d e+e~ 


^opt,N LO 


10 


0.08108 


0.08093 


20 


0.06574 


0.06565 


30 


0.05907 


0.05900 


40 


0.05508 


0.05503 


50 


0.05233 


0.05228 



Table 4: Same as in Table 3, but for n/ = 4. The PMS parameter in NLO is 
approximately r\ = —0.71. 

MS scheme and 0.05975 in the NNLO PMS. In Table 4 we give numerical values of 
the improved predictions in the PMS scheme, for several values of -y/i/A^-. We find 
that in the improved approach the NNLO PMS predictions are very close to NLO 
PMS predictions. We see therefore that also for rif = 4 the theoretical uncertainties 
in the improved predictions for 5 e + e - are very small. 

Finally let us consider the case of rif = 3. We have then = 5.23783. In 

(2) 

Fig. 13 and Fig. 14 we compare the plots of o e + e - as a function of r 1; for several 
values of C2, with a/s/A^^ = 9, obtained with the conventional NNLO expression 
(Fig. 13) and with the numerical evaluation of the contour integral (Fig. 14). Again, 
we find dramatic reduction in the RS dependence, despite rather low energy. It is 
interesting that in the rif = 3 case the addition of it 2 corrections given by Eq. (|35| ) 
and Eq. (|36|) does not result in the improvement of the conventional predictions. In 

(2) 

Fig. 15 we show the contour plot of 8 e + e - obtained from the improved expression for 
\fs/ A.jjs = 9- Similarly as for other numbers of flavors we obtain in the improved 
approach a very small variation in the predictions when parameters are changed over 

(2,) 

the I = 2 region of parameters appropriate for 5 D — the variation is of the order of 
0.8% of the PMS prediction 0.07756. (We have verified that this situation persists 
down to \/s/A|^ = 4.) The improved prediction in the MS scheme is 0.07719. 
For comparison, in the conventional approach we obtain 0.08097 in the NNLO PMS 
and 0.08244 in the NNLO MS scheme. In Table 5 we give numerical values of the 
improved predictions in the PMS scheme, for several values of y/s/A^^. With this 
results we conclude, that the rif = 3 NNLO expression for S e + e -, obtained by eval- 
uating the contour integral (|3^) numerically, has very small theoretical uncertainty, 
even for rather small values of \fsl A^—. This situation is similar to that found for 

v ' M §— 

the QCD corrections to the tau decay [29|, f|. 



(2) 

The behavior of o + at very low energies and the problem of existence of the 
fixed point in the improved approach would be discussed in a separate note 
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front, N N LO 
d e+e~ 


^opt,N LO 


5 


0.09624 


0.09421 


7 


0.08475 


0.08312 


9 


0.07756 


0.07619 


11 


0.07255 


0.07136 


13 


0.06879 


0.06774 



Table 5: Same as in Table 3, but for rif = 3. The PMS parameter in NLO is 
approximately r\ = —0.81. 

6 Summary and conclusions 

Summarizing, we have analyzed the RS dependence of the conventional NNLO ex- 
pression for 5 e + e ~ using a systematic method described in |T], 0, [J . We found rather 
large variation in the predictions. We have also investigated an improved way of 

(2) 

calculating o y + , which relies on a contour integral representation for this quantity 
and a numerical evaluation of the contour integral. We found that the stability of 
8 + - with respect to change of the RS is greatly improved when the contour integral 
approach is used. Also, in the improved approach the difference between optimized 
NNLO and NLO predictions was found to be much smaller than in the case of the 
conventional expansion. We conclude therefore that the theoretical uncertainties in 
the NNLO QCD predictions for S e + e - are very small, even at low energies, provided 
that large 7r 2 terms, arising from analytic continuation, are treated with due care. 
We observed that the optimized predictions for S e + e -, obtained in the contour inte- 
gral approach, lie in general below the predictions from the optimized conventional 
expansion. However, for nf = 5 the change in the fit of Aj^L to the experimental 
result came out to be small. 

Note added. After this paper was completed, a related work was brought to 



our attention ||50|1 , in which the RS dependence of the QCD corrections to the total 



hadronic width of the Z boson is discussed. In |hJ it is observed, that by using the 
contour integral to resumm the large 7r 2 contributions one reduces the scale depen- 
dence of the QCD predictions. This result is in agreement with our observations, 
since the dominant contribution to Y h £ d comes from expression identical to 5 



Let us note that the result reported in |5(J was anticipated already in |J. However 



the approach adopted in |5(| differs from our approach in several ways. The authors 
of []l)D| do not discuss the choice of the range of scheme parameters used in their 
analysis. In their investigation of the conventional expansion for r^ ad they use a 
smaller range of parameters than the one used above for n/ = 5. In particular, 



the PMS parameters are outside the range considered in |p0[ . In the analysis of 
improved predictions for Y h ^ d the authors of [50] limit themselves to the discussion 



of the renormalization scale dependence, fixing the /5-function to the MS value. 
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There is also a technical difference that authors of [|50| use approximate analytic 
expression for the running coupling constant to integrate along the contour in the 
complex energy plane, whereas we use exact numerical solution of the two or three 
loop renormalization group equation. 
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Figure Captions 



(2) 

Fig. 1 The contour plot of 8 e + e - as a function of the parameters r\ and c 2 , with 
rif — 5, for y/s/Aj^g = 75. The region of scheme parameters satisfying the condition 
( |TTD has been also indicated for I = 2 (the smaller region) and for 1 = 3. 

(2) 

Fig. 2 The maximal and minimal values of 8 + in the I = 2 (dash-dotted line) 
and I = 3 (dashed line) allowed regions, with rif = 5, as a function of y/s/Ajj^. 
The PMS prediction is also shown, and the experimental constraint 8 e ^? e -(\/s = 
31.6 GeV) = 0.0527 ± 0.0050 [55] is indicated for comparison. 



(2) 

Fig. 3 The contour plot of e + e - as a function of the parameters r\ and c 2 , with 
nf = 4, for y/s/Aj^g = 30. The region of scheme parameters satisfying the condition 
(|ITD with I = 2 has been also indicated. 

(2) 

Fig. 4 The variation in the predictions for e + e - when the scheme parameters are 
changed over the I = 2 region, with nf = 4, as a function of \/s/A^^. The up- 
per curve corresponds to r± = 3.10 and c 2 = 6.65, the lower curve corresponds to 
7*1 = —4.32 and c 2 = 0. For comparison the PMS prediction is shown. 

Fig. 5 Same as in Fig. 3 but for nj = 3 and y/s/ Aj^ = 9. 

(2) 

Fig. 6 The variation in the predictions for 8 e + e - when the scheme parameters are 
changed over the / = 2 region, with n/ = 3, as a function of y/s/A-^. The up- 
per curve corresponds to 7*i = 2.71 and c 2 = 5.71, the lower curve corresponds to 
ri = —3.21 and c 2 = 0. For comparison the PMS curve is shown. 

Fig. 7 5g+ e _ as a function of r 1; for several values of c 2 , forny = 5 and \/s/A^^ = 75, 
obtained with the conventional NNLO expression. For comparison also the NLO pre- 
dictions are indicated. 

Fig. 8 5g+ e - as a function of 7*1, for several values of c 2 , for n/ = 5 and y/s/Aj^ = 75, 
obtained with the numerical evaluation of the contour integral. For comparison also 
the NLO predictions are indicated, and the predictions obtained from the conven- 
tional expansion supplemented by the 0(a 4 ) and 0(a 5 ) corrections given by Eq. fl3"5p 
and Eq. ©. 

(2) 

Fig. 9 Contour plot of 8 e + e - obtained from the improved expression for rif = 5 and 
v/i/A^r = 75. The regions of scheme parameters satisfying the condition ( p4|) with 
/ = 2 (the smaller region) and I = 3 have been indicated, assuming p 2 = 

Fig. 10 Same as in Fig. 7, but for rif = 4 and y/s/A^L = 30. 
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Fig. 11 Same as in Fig. 8, but for rif — 4 and y/s/A^ = 30. 

Fig. 12 Same as in Fig. 9, but for rif = 4 and y^/A^ = 30. Only the I = 2 region 
has been indicated. 

Fig. 13 Same as in Fig. 7, but for rif = 3 and y/s/A^ = 9. 
Fig. 14 Same as in Fig. 8, but for rif = 3 and ^/s/A^ = 9. 

Fig. 15 Same as in Fig. 9, but for rif — 3 and t/s/Aj^ = 9. Only the I = 2 region 
has been indicated. 
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